Theory of Excitonic States in CaBg 
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We study the excitonic states in CaBg in terms of the Ginzburg-Landau theory. By minimizing 
the free energy and by comparing with experimental results, we identify two possible ground states 
with exciton condensation. They both break time-reversal and inversion symmetries. This leads 
to various magnetic and optical properties. As for magnetic properties, it is expected to be an 
antiferromagnet, and its spin structure is predicted. It will exhibit the magnetoelectric effect, and 
observed novel ferromagnetism in doped samples and in thin-film and powder samples can arise 
from this effect. Interesting optical phenomena such as the nonreciprocal optical effect and the 
second harmonic generation are predicted. Their measurement for CaBe will clarify whether exciton 
condensation occurs or not and which of the two states is realized. 



PACS numbers: 71.35.-y 75.70.Ee 77.84.Bw 

I. INTRODUCTION 

CaBg and its doped system Cai-ajLa^jBe have been 
the subject of intensive studies since the discovery of 
"high temperature ferromagnetism"EJ. The understand- 
ing of the parent compound CaBe, which we focus in 
this paper, is essential to that of this novel ferromag- 
netism. This ferromagnetism is peculiar in some respects: 
(i) the Curie temperature (~ BOOK) is high in spite of a 
small moment (0.07/LtB/La), (ii) it occurs only in a nar- 
row doping range (x < 0.01), and (iii) there are no par- 
tially filled d- or /-bands. Soon after this discovery, two 
scenarios are proposed; one is, based on a ferromagnetic 
phase of a dilute electron gaso, and another is based on 
an excitonic stateoa where the doped carriers are con- 
cluded to be fullvrnSpin-polarizedQ. Earlier band struc- 
ture calculationsoEl found the small overlap of the con- 
duction and valence bands near the three X-points. In 
the de Haas- van Alphen measurement, the Fermi surfaces 
of both electrons and holes are observedEl, which confirms 
the smallness of the gap. Furthermore, the symmetry of 
the wave function at the bottom of conduction band (the 
top of the valence band) forbids a finite dipole matrix ele- 
ments. Therefore the dielectric constant is not enhanced 
even when the band gap is reduceaJ, and the excitonic 
instability is not suppressed; CaBg can thus be regarded 
as a promising candidate for the excitonic state. A 
cent band calculation based on GW approximationll2 
on the other hand, predicts a large band gap of ~ 0.8 
eV, which is too large compared with the exciton bind- 
ing energy of the-prder of Tc ~ 600K. Measurements 
on photoemissiontJ and soft X-ray emissiorO observed a 
band gap of '--^ 1 eV, consistent with this GW calculation; 
yet a surface effect in photoemission should be seriously 
considered before any conclusion is drawn. Particularly, 
another GW calculatiorllj found the small overlap simi- 
lar to the LDA calculation. Furthermore, there are sev- 
eral experimental evidences showing a symmetry break- 
ing with magnetism below Tc 600K . X-ray diffraction 



and Raman scattering experimentJli shows an anomaly 
at Tc 600K, below which the symmetry is-Jowered from 
cubic to tetragonal. A /iSR measurementlla detected an 
internal magnetic field, which suggests an existence of a 
magnetic moment. 

These findings contradict the large band gap because 
any instabilities are unlikely if that is the case. In addi- 
tion, sensitivity to impurities and defects suggests a small 
band gap/overlap. Thus, it seems likely that the exci- 
ton condensation occurs in this compound due to a small 
band gap/overlap. One can then attribute the high Curie 
temperature in Cai-ajLa^jBg to the excitonic instability 
in the parent compound, which is determined by the ex- 
citon binding energy. Nonetheless, even if we assume 
the exciton condensation, there remain many mysteries. 
The peculiarity of ferromagnetism in Cai-ajLa-ajBg reveals 
itself in the strong sample dependence. In Ref. ^ the 
doped carriers are assumed to be essential to ferromag- 
netism, but recent experiments could not find the corre- 
lation between the transport properties and magnetism. 
The smallness of the magnetic moment (~ O.OT/is/La) 
is not fully accoup±jed,iQpr| although numerous attempts 
have-been madedllj'EIlHE^. Moreover deficiency iruCa 
sitjescS or doping of divalent elements such as Ba£J or 
SrEj also induces ferromagnetism. Another hint is that 
the fermmagnetic moment is enhanced in the thin-film 
sampleE3, or -near the surface from electron spin reso- 
nance (ESR)lj. This fact suggests that some kind of 
symmetry lowering is related to the ferromagnetism. 

Because of this controversial situation, it is worthwhile 
to give solid theoretical results which are independent of 
the details of microscopic models. In this paper we deal 
with this problem from the viewpoint of symmetry. The 
only assumption we take as the basis of our analysis is 
that the parent compound CaBg is a triplet excitonic 
insulator with the order parameters made from the Xg- 
(conduction band) and X3- (valence band) wavefunctions 
at X-points. The analysis based on the Ginzburg-Landau 
expansion becomes exact by considering the symmetry 
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properties of the order parameters. All possible states 
are classified according to the irreducible representation 
of the magnetic point group. We can then predict many 
nontrivial physical properties for each of the states, which 
can be tested experimentally. Most important of these 
are that these states break both the time reversal (R) 
and space inversion (/) symmetries, while their prod- 
uct RI remains intact. This implies that CaBg is an 
antiferromagnet (AF). This also leads to the magneto- 
electric (ME) effect; namely the ferromagnetic moment 
is induced when the electric field is applied or the electric 
polarization is induced when the magnetic field is applied. 
This offers a new mechanism of the novel ferromagnetism 
in Cai-a^LarcBe, i.e. ferromagnetism is induced near im- 
purities or surfaces by the ME effect We note that the 
present theory of ferromagnetism is different from those 
in Ref. |[^, though all these theories assume exciton con- 
densation. The relation between our theory and the the- 
ories in Ref. is presented in det;ui in the last section. 
We have published a short versionEa of this paper, but 
this paper contains new results such as a pattern of the 
magnetic moments, and optical properties, as well as full 
and detailed description of the discussion. The plan of 
this paper is as follows. In Section II, the excitonic order 
parameters are defined and the Ginzburg-Landau expan- 
sion is developed in III. Physical properties of each state 
is described in Section IV, and V is devoted to discussions 
and conclusions. 

Here we remark a role of spin-orbit coupling. In sub- 
sequent discussions we take into account the spin-orbit 
coupling. Its presence is crucial in our analysis by the 
GL theory. Limit of zero spin-orbit coupling is discussed 
in Appendix |d|, and it can be contrasted with the nonzero 
case. In our numerical estimation for coefficients in GL 
free energy (Fig. ||), we neglected the spin-orbit coupling, 
because it greatly simplifies the calculation; in interpret- 
ing its results, however, we should keep in mind that we 
are always treating the case of nonzero spin-orbit cou- 
pling. 



II. EXCITONIC ORDER PARAMETERS 

Let bka and a^^ denote annihilation operators of elec- 
trons with snin a at the conduction and the valence bands 
respectivelyo. When excitons are formed, excitonic order 
parameters {b^^ak'fj) will have nonzero values. ExcitoHic 
instability occurs only in the vicinity of the X pointsBLl, 
i.e. = (7r,0,0), Qy = (0,^,0), = (0,0,^) in 

the cubic Brillouin zone, where the conduction and the 
valence bands approach each other with a small over- 
lap/gap. We make the following assumptions; 

(i) The order parameters are fc-independent near the 
X points. 

(ii) The order parameters connecting different X points 
are neglecteda. 



Consequently, we keeft only the order parameters 
v{Q^)a|3 = (5^^„aQ./5)§. Still wc have 2 • 2 • 3 = 12 
complex order parameters involved and the problem is 
quite complicated. The GL theory is so powerful that it 
helps us from this difficulty. The cubic (Oh) symmetry 
of the lattice considerably restricts the form of the GL 
free energy <&, as is similar .in . the GL theory for uncon- 
ventional superconductorfcSEZI. In this section and in the 
next section, we construct the GL free energy and discuss 
its mimima. The argument is rather technical, because 
we will make full use of the point-group symmetry to deal 
with the twelve complex order parameters. Readers who 
are not accustomed to detailed point-group discussions 
can jump to Section IV. 

At each X point, the fc-group, i.e. the group which 
keeps k unchanged, is tetragonal (D4h), and the con- 
duction and the valence band states belong to Xg 
and X3 repjjfisentations, respectively according to band 
calculationBCI, when we take the origin as the center of 
a Bg tetrahedron. In the absence of the spin-orbit cou- 
pling, the order parameter 'q{Qi)ap transform as a 
irreducible representation X3 x X3 = X'j^. From now on 
we shall take the spin-orbit coupling into account. Then, 
the representation of 77 will be altered. We restrict our 
analysis on the triplet channel for the excitons, because 
the exchange interaction usually, favors the triplet state 
compared with the singlet stateB. Then, the spin-1 rep- 
resentation is multiplied to X'j^ , and the triplet order pa- 
rameters follow the representations X4 -I- Xg in the D4h 
group. Let us now construct the order parameters ex- 
plicitly. We take the Qz point as an example; the triplet 
order parameters have three components, following the 
spin-1 representation. By recalling that a spherical har- 
monics with ^ = 1 is represented by x ± iy, z, we can 
recombine the triplet order parameters as 



VyiQz) 

VziQz) 

or in a short form 



V2 
1 

'V2 



{miiQz) + vi^iQz)), 



a,/3 



(1) 

(2) 
(3) 

(4) 



The reason why we introduced these linear combina- 
tions is because they are convenient for subsequent sym- 
metry discussions. Their phase factors are chosen so 
that the time-reversal R operates as complex conjugation 
Ri]i{Qa) — ViiQa)- By investigating their transforma- 
tion property under D4h, which transforms Qz onto itself, 
we find that r]z follows the X'^ representation, while rj^ 
and r]y follow the X'^ in D4h. Here the spin-quantization 
axis in (P-(||) is taken to be the z-axis. Since we are tak- 
ing into account the spin-orbit coupling, these uo- and 
down-spins should be interpreted as pseudospinscZI. 
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When we consider the three X points, the cubic sym- 
metry is restored. We find that nine components of the 
order parameters are classified into irreducible represen- 
tations of Oh as Fis -I- Fis -I- r25, all of which are 3- 
dimcnsional. Let us call basis functions as 77i(ri5,l), 
77i(ri5,2) and ?7j(r25), where i = x,y,z and j = 
x(y^ — z'^),y{z'^ — x^),z{x^ — y^). These basis func- 
tions transform like a vector (a;, y, z) for Fis and like 
(a;(y^ — z'^),y{z'^ — x'^),z{x^ — y"^)) for r25. The basis 
functions are given as 



- \{KvniQy) + vii{Qy)) + {vniQ^)-vi'\{Qx))} 
\{KmiiQy) + vi]{Qy)) ~ {vniQx)-vuiQx))}- 

We defined these linear combinations in order to facilitate 
the construction of the GL free energy $ in the next 
section. Here, the spin-quantization axis in r]af3{Qi) is 
taken as the -|-i-axis (z = x,y,z). Other components are 
obtained by cyclic permutations of x,y,z. It is easily 
shown that Irj ~ rj, and Rr] = rj* , where we defined 
for brevity 

r7(ri5, 1) = (77.(ri5, 1),r7y(ri5, l),Vz(Ti5, 1)), 
r7(ri5,2) = (77,(ri5,2),77^(ri5,2),7;,(ri5,2)), 

»7(r25) = {Vx{y^-z^)i^25) I Vy(z'^-x^) (r25), 

?7z(a;2-y2)(r25)). 



III. GINZBURG-LANDAU THEORY AND 
DETERMINATION OF THE GROUND STATE 



A. Quadratic Order Terms 

Let us now write down the GL free energy in terms of 
these order parameters. The GL free energy $ should be 
invariant under the elements of Oh and under the time- 
reversal R. We make two remarks helpful in writing down 
First, only even-order terms in rj are allowed by the 
inversion symmetry. Second, owing to the time-reversal 
symmetry, the order of Imr/ in each term should be even. 



Thus, $ is given up to quadratic order as 
$(2) =Ai^(Re,7,(ri5,l))' 

i 

+ A2^(Rer;,(ri5,2))2 

i 

-f Aa^Re 77,(ri5,l)Re r;,(ri5,2) 

i 

+ A4^(Rer,,(r25))2 

j 

+ Si ;^(Im 77,(ri5, l) f + B2 5](Im //.(Lis, 2))^ 

i i 

+ S3 ^Im 77, (Pis, l)Im ?7,(ri5, 2) 

i 

-^B4^(Imr;,(^25))^ (5) 
i 

where Ai , Bi are constants. As there exist no other sym- 
metries than considered above, these Ai and Bi are in- 
dependent parameters. 

Let us determine the ground state of the system from 
the GL free energy. Minimizing we find that one of 
the following states will be realized as the temperature is 
lowered; 

(A) Rcr7(ri5,l) = cRer7(ri5,2)^0, 

(B) Re r7(r25) ^ 0, 

(C) Im77(ri5,l) = clmr7(ri5,2)^0, 

(D) Im r7(r25) ^ 0, 

where c is a nonzero constant. Therefore, condensations 
of excitons in Lis and in r25 do not occur simultaneously. 

In all cases (A)-(D), all directions of the vector ij are 
degenerate, and this degeneracy will be lifted in the quar- 
tic order terms in as we shall see afterwards. All these 
states are accompanied by a lattice distortion. This dis- 
tortion, however, is expected to be rather small, because 
it couples to the order parameter in the quadratic order, 
not linear order. Indeed, the distortion detected by tiie 
X-ray scattering is tetragonal and is as small as 0.03%O. 
This smallness of distortion is a common feature also in 
other hexaborides; for example, there observed no lattice 
distortion associated with low-temperature antiferromag- 
nctisip-|in GdBgES, and antifcr|Mquadrupolar ordering in 
CeBg". Ferromagnetic EuBgE^J should also have a dis- 
tortion away from a cubic lattice, though not observed 
experimentally. It is because ferromagnetism cannot oc.- 
cur in a cubic lattice from general symmetry argumentcil. 

B. Microscopic Calculation for Quadratic Order 
Terms 

To study which state is realized, we shall calculate tlie 
coefficients in by the Hartree-Fock approximationB, 
including coupling terms between the order parameters 



4 



corresponding to the different X points. Because the 
spin-orbit couphng is weak in this compound, we shaU 
neglect it in the following microscopic calculation. 

In order to derive the quadratic term of the GL 
theory for the triplet order pararoeters, we have only to 
concentrate on the following tcrmLl. 



boron sites, as shown in Fig. |i| (C). As a simple estima- 
tion, we approximate the orbitals of the boron and the 
calcium ions as those of a hydrogen-like atom with charge 
+2e and +3e, respectively. We believe that using more 
realistic orbitals in the solid will not change our conclu- 
sions described below. Owing to the time-reversal sym- 
metry, the wavefunctions at each X-point can be taken as 
real. Still, there remains uncertainty in relative signs of 
the Bloch wavefunctions among the X-points. We define 
this as follows; 



where 



V, 



da; J da;>^^ (a;)^^ (x) 0^" {x'Wq[ {x'), 



and e is a dielectric constant, 0q. (a;) is the Bloch wave- 
function of the K-band {k = a, 6) at fc = Qi, taken as 
real. We assume that screening of Coulomb interaction 
is weak. Let us construct the wavefunctions involved in 




FIG. 1: (A) Wavefunction of X3 from the d-orbit of the Ca- 
(B) Wavefunction of X3 from the p-orbits 



sites, 



iCa:d — orbit 



of the B-sites, <j)^'^~f^''' (C) Wavefunction of X3 from the 

p-orbits of the B-sites, '^q^x™'"*- wavefunctions de- 

picted in this figure have crystal momentum Qz = (0, 0, n) 

the exciton formation. The wavefunction for the conduc- 
tion band (j>Q, (x) at the X points mainly consists of the 
p-orbits of the boron ions and the d-orbits of the cal- 
cium ions. These constituent wavefunctions, </'g''x"'^'* 
and <^q'J;17°'''", are depicted in Fig. (A) (B). The wave- 
function (j)Q, (x) is a bonding orbital of these two: 



B:p — orbit _|_ — j,^0Ca:d — orbit 



(7) 



On the other hand, the wavefunction of the valence band 
at the X points is mainly composed of the p-orbits of the 



where Cg^^^' is a three-fold rotation around the [111] di- 
rection. This particular choice of relative signs of the 
wavefunctions at different X-points does not affect phys- 
ical consequences given below; another choice leads to 
the same physical consequences. 

Now we can evaluate the coefficient of eq.(p|) from 
eq. (||) . Because there are some relation such as Vq°^q^ = 

^Qyh. " ^Qyh^ = ^qf Qx ' ^6 have only to evaluate 
three quantities a = (3 = V^l%^, 7 = V'^^t'^Q^- 

In terms of these quantities, the coefficients in the GL 
free energy (H) are expressed as 



Ai^C + a, A2^C + a + l3-j, 

A3 = 2V2(/3-7), Ai^C + a~l3 + j, 

Bi = C-a, B2^C-a-P-j, 

B3 = -2V2(/3 + 7), S4 = C - a + /3 + 7, 



(8) 



where C is a constant. These satisfy relations (Dl) in 
Appendix which reflects the SU(2) symmetry in the 
spin space in the limit of zero spin-orbit coupling. To 
determine which type of excitons will condense, we diag- 
onalize as 



$(2) = A5°|Rer,(Fi5,+)P + Ai7|Imr,(Fi5,+)P 
+A^'=|Rer7(ri5, + Ai™|Imr,(ri5, -)\' 



(9) 



where 



= 



Air = 



i,Re 
\ Im 

^r2t 



C-f -(2a + /3-7±3(/3-7)), 
C+i(-2a-/3-7±3(-/3-7)), 

C - a + /3 + 7 = A^!", 



(10) 



and 



77(ri5, +) = ^(»7(ri5, 1) + \/2r,(Fi5, 2)), (11) 
r7(ri5, -) = -^(V2r,(Fi5, 1) - »7(ri5, 2)). (12) 
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The degeneracy in Ap™ = A'f'^ Ap°^ — is a conse- 
quence of (Dl) which is vahd in the absence of the spin- 
orbit couphng; this degeneracy is Hfted by the spin-orbit 
couphng. 

Let us turn to numerical evaluation of a, (3, 7. a, /3, 7 
can be simplified as follows, 



unit cell 



iGy 



1 



Go + G|2 



unit cell 



unit cell 



1 



i(GQ+G)-y 



unit cell 



unit cell 

where Go — *(7r, tt, 0), and the summation over G rep- 
resents the one over reciprocal vectors. The function 
4'Qi (^) is a localized atomic wavefunction, and we ap- 
proximate it as 



N<I)q,[x) : inside the unit cell, 
: outside the unit cell. 



(13) 



where T denotes a translational vector of the lattice 



Hence, it satisfies (j)Q,{x) 



^(a;-r)e'"^-Q% 
We 

have omitted the factor e^/e in the expressions of a, /?, 7, 
since it is irrelevant for the present discussion. Shown in 
Fig.^ are the coefficients of eq. (^) as a function of c de- 
fined in (^. The constant part C in (10) common among 
these coefficients is omitted in this calculation, since it 
is unnecessary for comparison of A^™ and A^° . As is 
shown in Fig.0, the coefficients involving the imaginary 
parts are smaller than those for the real parts. Therefore, 
the states (C) and (D) would be favorable than (A) or 

(B) . In the notation in Ref. ^ the type-II and type-Ill is 
more favorable than the type-IV and type-I. The states 

(C) (D) breaks the time-reversal symmetry. Thus they 
are magnetic states, which seems to fit an appearance 
of ferromagnetism in several conditions like powder or 
thin film in CaBg. On the other hand, the states (A)(B) 
preserve the time-reversal symmetry and thus nonmag- 
netic. From symmetry consideration, this implies that 
neither the ME nor the piezomagnetic (PM) effect will be 
observedE3. Roughly speaking these states are far from 
showing ferromagnetisrn—. Thus, hereafter we shall con- 
centrate on (C) and (D)cj. 

Figure || shows A^f^ > = Ap™ , meaning that the 
condensation of excitons in Fis and those in seem 
to occur simultaneously. This is an artifact of approxi- 
mation of zero spin-orbit coupling. We note that we are 
working with the presence of spin-orbit coupling; these 
two condensations do not occur simultaneously, and ei- 
ther (C) or (D) is selected. 
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FIG. 2: Coefficients Xf", \^]^ of the GL free energy for real 
and imaginary parts of the order parameters, a =4. 15 A is 
the lattice constant of GaBg. The parameter c controls the 
X3 wavefunctin of the conduction band. It is defined in (^) 
as a coefficient of mixing of the p-orbit of the boron and the 
d-orbit of the calcium. 



C. Quartic Order Terms 

Let us consider quartic order terms in the GL free 
energy, which lifts the degeneracy in the direction of rj. 
The terms in containing only the imaginary parts of 
the order parameters arc written as 

$(4) = 

Dl ^^(^15, l))^! + El ^(Im 7/,(ri5, 1))^ 
+ D2 77»(ri5, 2))2| + ^2 ^(Im 77,(ri5, 2))4 

+ D3 ^.(r25))'| + ^3^(Im r7,(r25))" 

+ Fi |E(I"^ '?^(ri5, i))^! |E(i™ ^'(^15, 2))2| 

+ F2 {3(Im 7/,(ri5, 1)2 - Im r;y(ri5, 1)^) 

X (Im r?,(ri5,2)2-Im 7jy{ri5,2f) 
+ (2Im 77,(ri5, If - Im r;,(ri5, 1)' - Im 77^(ri5, 1)') 

x(2Im r7,(ri5, 2f - Im 77,(ri5, 2)^ - Im 7^y{ri5-2f)} 
+ F3 

X (Im 7]x{Ti5, l)Im r7y(ri5, l)Im ?7^(ri5,2)Im 77j^(ri5,2) 
+ Im riy{Ti5, l)Im %(ri5, l)Im 77j,(ri5, 2)Im 77^(ri5, 2) 
+ Im 77^(ri5, l)Im rix{Ti5, l)Im 77^(ri5,2)Im T]a:(Ti5,2)) 

+ ^Imr;,,(ri5,l)Im77,(ri5,2) 
i 

X I Glim r/,(ri5, 1)2 + G2 ^Im 77,,(ri5, 1)^ 
+G3lm r;, (Fis , 2)2 + G4 ^ Im r;,, (Fis , 2)^ | . 
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Here coupling terms between Fis and are not writ- 
ten, because we are dealing only with condensation of 
one type of excitons. It is noted here that the coupling 
of the order parameter rj and the lattice distortion u such 
as u^T]^ also produces the quartic terms in rj after inte- 
grating over u, and hence lifts the degeneracy. This effect 
is included in By minimizing <i>(^) -|- we find 

four possibilities; 



(C-1) Im J7(ri5, 1) = cilm v{^i5,2) = (0, 0, C2), 

(C-2) Im J7(ri5, 1) = cilm ^(^15, 2) = (ca, C2, ca), 

(D-1) ImJ7(r25) = (0,0,C2), 

(D-2) Im J7(r25) - (C2,C2,C2), 



where c's are constants. There are three equivalent direc- 
tions [001] for (C-1) and (D-1) and four directions [111] 
for (C-2) and (D-2); we take ony one of them without 
loss of generality. The direction of the lattice distortion 
is tetragonal in (C-1) and (D-1) and is trigonal in (C-2) 
and (D-2). We remark here on thexecent results of X-ray 
scattering and Raman scatteringliJ, which strongly sup- 
port our theory. They show a tetragonal distortion below 
600K, which indicates that (C-1) or (D-1) is realized in 
CaBg. Furthermore, inversion-symmetry bijeaking does 
not manifest itself in the Raman spectrumEfl, which is 
consistent with (C) or (D), but not with (A) or (B). The 
reason is the following. In all cases (A)-(D) the inver- 
sion symmetry is broken. In (A) and (B) this broken 
inversion symmetry will appear in the Raman spectrum. 
The cases (C)(D), on the other hand, are invariant under 
an RI symmetry, as we shall see later. Since the phonon 
spectrum is time-reversal invariant, it is automatically in- 
variant under RI - R = I, i.e. inversion symmetric. Thus, 
although (C)(D) breaks inversion symmetry, it does not 
appear in the Raman spectrum. Therefore, the phonon 
spectrum, which is even in time-reversal, automatically 
become even in inversion. Thus, although (C)(D) breaks 
inversion symmetry, it does not appear in the Raman 
spectrum. 

To summarize, (C-1) and (D-1) are the only possi- 
bilities totally consistent with the X-ray and Raman 
scattering^. In the states (C-1) and (D-1), the con- 
densed excitons in each case are given as 



where Ci are real. To be explicit, we can rewrite them as 
(C-1) : (b^Q^^aQ^i) = -{b^Q^^aQ^^) = ^, (16) 

= ^{bQ,]aQ^i) = -i(foQ^|aQ,T> = (17) 
(D-1) : {bQ^i^aQ^l) = (&Q„iaQ„T> 

t ' ^6 

= iiVQ^^aq^i) = -i(6^j^|aQ,T) = (18) 

These expressions are too complicated to extract in- 
formation for physical properties; instead our approach 
is based on symmetry, namely, magnetic point groups. 
Magnetic point groups G for these states are easily ob- 
tained, if we consider which symmetry operations in- 
cluding R keeps the values of the order parameters un- 
changed. The results are 

(C-1): G = 4/m'mm = C4V x {E, RI} 

— {E, C2Z adb} X {E,RI}, 

(D-1): G = 47m'm'm = D2d x {E, RI} 

~ {E, Sf^, C2z,G2a,C2b, CFx,(Jy} X {E, RI}. 

Here symmetry operations are defined as follows; E: 
identity; C2j (j = x,y,z); two-fold rotation around the 
j-axis; C2a, C2b- two-fold rotation around the [110] and 
[110] axes, respectively; C|^: ±7r/2 rotation around the 
z-axis; aj {j = x,y,z): reflection with respect to the 
plane normal to the j-axis; S^^ = Cf^I; <7da, o-db- re- 
flection with respect to the (110) plane and the (110) 
plane, respectively. These will be used in the next sec- 
tion to make various predictions of the present material. 
In (C-1) and (D-1), when we fix the axis of tetragonal dis- 
tortion, there are two types of degenerate AF domains, 
related with each other by the time-reversal. Since there 
are three choices of axes for the tetragonal axis, total de- 
generacy is six in (C-1) (D-1), which is equal to an order 
of the quotient group (Oh x {E,R})/G. We can draw 
some analogies with anisotropic superconductivity (SC). 
The order parameters t] of triplet excitons correspond to 
the d-vector in triplet SC. It is nevertheless misleading 
to look for SC counterparts of our phases (C-1) (D-1), 
because our order parameters are confined in the neigh- 
borhood of the three X points. They are triplet and even 
functions in k, which never occurs in the SC. 



IV. PREDICTION OF PHYSICAL PROPERTIES 



A. Magnetic Properties 



(C-1) : r,(Q,) = J7(Q^) = (0,0,ic4), 

r,(Q,)- (0,0, ics), (14) 

(D-1) : viQ.) = -ViQy) = iO,0,ice), 

V{Qz)^0, (15) 



A crucial observation for prediction of physical prop- 
erties of these states is that they do not ba-eak the RI 
symmetry, while both R and / are brokenu. There are 
a number of compounds, such as CraOa, known to pos- 
sess this feature, leading to several magnetic properties 
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listed below. First, the RI symmetry prohibits a pres- 
ence of any uniform magnetization. Dzyaloshinskii used 
this symmetry to explain why weak-ferromagnetism is 
present in Q;-Fe203 while not in Cr203E3. Thus the states 
(C-1) and (D-1) are antiferromagneiic. This agrees with 
the result of the fiSK measurementEj with a moment of 
0.0039/XB/mol. Note that the magnetic unit cell is iden- 
tical with the original unit cell. Thus, no extra Bragg 
spots appear below the AF phase transition. Second, 
the i?I-aymmetry also prohibits the piezomagnetic (PM) 
effect a uniform stress cannot break the RI symmetry 
which hinders ferromagnetism. Nevertheless, a gradient 
of stress can break this symmetry and will induce ferro- 
magnetism. 

Third, the i?/-invariance rcsiilta-jn the linear ME ef- 
fect, as is first observed in Cr203EZH23. The i?/-invariance 
allows the free energy to have a term like GijHiEj, caus- 
ing the magnetization M proportional to E, and the 
polarization P proportional to H: 

M = GE, P = ^GH. (19) 

Roughly speaking, this occurs because an external elec- 
tric field E breaks this RI symmetry and enables 
ferromagnetismO. In the GL language, this can be stated 
as follows. An electric field E belongs to the Fis repre- 
sentation, and couples linearly with the order parameters 
in (C-1) as 

5$ = -C\E ■ Re r7(ri5, 1) - C2E ■ Re J7(ri5, 2), (20) 

in the lowest order. Here the imaginary parts of the or- 
der parameters are absent due to invariance of (5<f> under 
time-reversal. Thus, in the presence of E, both the real 
and imaginary parts of the order parameters acquire non- 
vanishing values; this breaks the RI symmetry, resulting 
in a ferromagnetic moment. As for the (D-1), similar 
effect can be found in 

3 

= -Ci ^(Im77,(r25)'(3Re77,(F25)K, - Re»7(F25) • E)) 

i=l 

-C2(Im77i(r25)Im772(F25)(ReJ7(r25) X E)3 

+Imr?2(F25)Im7;3(r25)(Rer7(r25) x E)i 
+Imr?3(F25)Im7;i(r25)(Rer7(r25) x E)2), (21) 

where T7(r25) — (?7i, ?72, %)■ Below the exciton condensa- 
tion temperature, Imr7(F25) has a nonvanishing value, 
which in turn brings about a linear coupling between 
Re?7(F25) and E. 

These properties described above are deduced solely 
from the i?/-invariance, and are common in (C-1) and 
(D-1). Meanwhile, detailed magnetic properties vary 
among them. Let us first investigate a AF magnetic 
structure within the unit cell. We can calculate pos- 
sible magnetic structure consistent with each magnetic 
point group. We cannot determine the quantitative dis- 
tribution of magnetic moments. Their distribution can 




FIG. 3: Schematic picture of spin structure for (C-1) and (D- 
1). The right panels in (C-1) and (D-1) show the top views, 
i.e. the views from the [001]-direction. 

be novel, since in the isostructural compound CeBg, po- 
larized neutron scattering shows that magnetic moments 
are at three locations: at the Ce sites, at the centers of the 
triangular plaquettes of Bg octahedra, and at the pidst of 
the B-B links connecting neighboring octahedra.E3 Any- 
way, to look at the difference intuitively, we displayed in 
Fig. 1^ directions of magnetic moments both at the boron 
sites and at the centers of the triangular plaquettes on 
Bg octahedra with the magnitudes arbitrary chosen. Let 
us call the six moments at each boron site as S±x, S±y, 
S±z', the obtained magnetic structure is 

(C-1) : -±(0,ai,0), S±j, = T(ai,0,0), 

S±, = 0, (22) 

(D-1) : = ±(0,ai,0), ^±^ = ±(01,0,0), 

S±, = 0, (23) 

where oi is a constant. On the other hand. Let SaP'j 
(a, /3, 7 = ±) denote a moment on a triangular plaquette 
with vertices Sax, S/3y, and S^z- Then they are given by 

(C-1) : S±±± = ±a3(-l, 1,0), S^±± = ±a3(-l, -1, 0), 
5±T± = ±03(1, 1,0), S±±^ - ±a3(-l, 1, 0), 

(D-1) : S±±± = ±(03, 03, -02), S^±± = ±(a3, -03, 02), 
S±^± = ±(-a3, 03, 02), S±±^ = ±(03, a3, 02), 

where are constants. These can be distinguished from 
each other by neutron scattering experiments. 

The difference in magnetic point groups is also re- 
flected in difference in a possible ME effect. The ME 
property tensor G in the term GijHiEj in the free enoMy 
is easily written down for each magnetic point groupE3: 



/ Gi 
(C-1) -Gi 
\ 
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in the original Cartesian coordinates of the cubic lattice. 
The property tensor G is antisymmetric in (C-1), which 
implies M±E and PJ-B. Thus, measurement of the 
ME effect with a single crystal of CaBg will reveal which 
of these two cases is realized. Experimentally, there are 
two types of AF domains, and the sign of the property 
tensor G is reversed when the staggered magnetization 
is reversed. In measurement of the ME effect, alignment 
of the domain structure is necessary. It is accomplished 
by means of a magnetoelectric annealing, in which the 
sample is cooled under both electric and magnetic fields. 
Experimentally, CaBg is not insulating, though the sam- 
ple dependence is rather large. Therefore, we cannot ap- 
ply electric field into the sample; we can neither perform 
magnetoelectric cooling nor the measurement of magne- 
toelectric effect. We should then resort to the optical 
measurements described in the next section. 

Domain boundaries between the two AF domains can 
exhibit interesting properties. As in boundaries between 
two SC domains with broken time-reversal symmetrytZl, 
localized current and magnetic moment are induced near 
the boundary. In the present case it is interpreted as the 
ME effect. 

It is reported_that weak diamagnetism is often ob- 
served in CaBgcrEj. It can be attributed to orbital mo- 
tion, as in the Landau diamagnetism. In semimetals like 
Bi and narrow-gap semiconductors, such diamagnetism 
appears. In Bi in particular, it is proposed as orbital 
diamagnetismc3. Hence, it is no wonder that the small 
band gap/overlap in CaBg originates diamagnetism. In 
reality, it is easily hidden by the AF or by the ferromag- 
netism near the surface. 



B. Optical Properties 

Optical measurements are known as powerful tools in 
identifying the magnetic properties of the material. In 
optical measurements, there are wide range of choices 
for tuning parameters, such as polarization and incident 
direction of light, providing us with lots of information 
of the material. One of the most spectacular examples 
is Cr203, where the second-harmonic generation (SHG) 
is utilized to get the image of the antiferromagnetic do- 
mains for the first timec2l. As shown in this example, 
magnetoelectrics such as Cr203 generally exhibit novel 
optical properties. Since our theory predicts that the 
present material CaBg is among the magnetoelectrics, 
its magnetic symmetry can be checked by various optical 
experirngnts . For example, nonreciprocal (NR) optical 
effectsOLJa occur in the magnetoelectrics, since tlie,ME 
tensor appears in the formulation of the NR effectalfl. In 
other words the ME effect corresponds to a static limit 
of the optical NR effects. Possible NR effects for all the 
magnetic point groups of magnetoelectrics are calculated 
in Ref. ^ for transmission and in Ref. ^ for reflection. 
Let us apply the essential results in these references. Here 
we note a difference in notation from Ref. for (D- 



1) 4'/m'm'm the x and y axes are rotated by 45 degrees 
around the z-axis from those in Ref. 

1. Optical Nonreciprocal Effects 

According to Ref.^, the multipole expansion of elec- 
tric and magnetic fields gives rise to various tensors, each 
of which is responsible for various optical effects. To 
the order of electric quadrupoles and magnetic dipoles, 
nonvanishing tensors in an _R/-invariant system such as 
CaBg are the polar i-tensor aQ/3(a;), the polar c-tensor 
"q/37('^)(= ^'a-ypi^))^ ^^'^ ^hc axial c-tensor Gafsito) in 
the notation of Ref. where aa/sij-^) is the polarizabil- 
ity tensor, while a^^^(a;) and Ga/3(w) are assigned to the 
NR optical effects. Note that the static limit Gap{i^ = 0) 
is just the ME property tensor Ga^- Magnetic symme- 
try determines nonvanishing components of these ten- 
sors. We can then calculate the refractive index using 
these tensors. The detailed derivation is developed in 
Appendix ^ and here we present only the results. Let 
n — n' + in" denotes the complex refractive index, where 
the real part n' is the refractive index and the imaginary 
part n" is the absorption coefficient. We also define a as 
(T = k/k, i.e. a normal vector along the incident direc- 
tion of light. The refractive index n for light along the 
±z-direction cr = (0,0, CTz) = (0,0, ±1) is obtained for 
two types of polarization as 

(C-1) : 

[E\\x, E\\y] -n^^ (^(JzfioCs + ^J ^JilCl + Ajiot^ , 
(D-1) : [E\\x] ■■ n=^ (a.fioC^ + filCl + ifioe^'^ , 

[E\\y] ■■ n='- (-UziiQC^ y^MoCf + Voe^^ 

where C5 — 2Gxy+ioa'^^^. Note that since Gafi and a'^^^ 
are c-tensors, C5 changes sign for different AF domains, 
which are related with each other by the time-reversal. 
Thus, for a single domain, (C-1) shows a directional bire- 
fringence, but there is no anisotropy in polarization, i.e. 
no birefringence in the usual sense. In (D-1), on the 
other hand, directional birefringence and usual birefrin- 
gence both exist. Hence, linearly-polarized light can dis- 
tinguish between these two cases. In fact we do not even 
need to fix polarization; unpolarized light is sufficient to 
distinguish between these two possibilities. Unpolarized 
light, which is an incoherent superposition of both x- and 
y-polarized lights, shows a refractive index of an average 
of two refractive indices for two types of polarizations 

(C-1) : n - ^ (^-a,mC5 + ^MoC|+4Moe.) , (24) 

(D-1) : n^yf4Ci+ifioe.. (25) 

Therefore, only in (C-1) does the system show directional 
birefringence for unpolarized light along the tetragonal 
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axis. This effect, i.e. directional birefringence-.for unpo- 
larized lights, is called a magnetochiral effected. This ef- 
fect is more prominent in the high-frequency region as esj- 
idenced in the X-ray measurement of Cr-dopcd V2O3EZI. 

Experimentally, domain structure is formed, which 
might obscure the experimental results. If we assume 
that the sample is a single crystal with perfect cubic 
structure above the Neel temperature T/v, there are two 
sources of domain formation below T^r. One is (a) the 
AF domains of opposite staggered moments, while the 
tetragonal axis are common. The other is (b) the do- 
mains due to different direction of tetragonal distortion 
below Tjv. This domain formation will affect the NR 
optical effects in the following way. For the (a)-type do- 
mains, the parameters C5 has different signs for different 
AF domains, which smears out the NR effects. Thus, the 
thin-film sample should have a single domain in the di- 
rection of light propagation. On the other hand, for the 
(b)-type domains, we can investigate the problem by set- 
ting the incident direction as x ot y direction in the above 
equations. The result for ex — {a^, 0, 0) = (±1, 0, 0) is 

[E\\y] ■■ n=^Je^h^, 

[E\\z]: n=y^(^l + ^C7) \ 

for both (C-1) and (D-1), where we set C7 = —Gxy + 
^ijj{a'zxx + ^'xxz)- 1^0^' the incident direction parallel to 
the y-direction, the result is similar, with E\\y and E\\x 
being interchanged in the above result. These imply that 
for cr 1 1 a; and <j\\y the NR optical effects do not take place. 
Thus only the light traveling parallel to the z-axis under- 
goes the NR optical effect; hence the (b)-type domains 
will not hinder the identification of ground states by the 
NR effect, although it is reduced by the factor of 3. 

In reflection, on the other hand, within the method 
in Ref. the reflection matrix R for light propagating 
along the z-axis defined by 

{E,)j = Rjk{Ei)k (26) 

is proportional to identity for both (C-1) and (D-1), 
where Ei and E^ are the electric fields for incident and 
reflected lights, respectively. This implies that there is 
no optical NR effect in reflection normal to the xy-plane. 
Nevertheless, since this effect is not prohibited by sym- 
metry, this NR effect in reflection can emerge in higher 
order than in Ref. ^ We note that in reflection ex- 
periments, the surface inevitably affects the spectrurrO. 
Thus, even if the NR effect in reflection are observed, it 
might be difficult to separate the result into the surface 
and bulk contributions, as it is controversial in Cr203. 

2. Second Harmonic Generation 

Second harmonic generation (SHG) is recently emerg- 
ing as a new and powerful tool for determination of mag- 
netic symmetry of the material. This effect appears both 



in transmission and in reflection. For separating bulk 
contribution from surface one, transmission measurement 
is preferable. At the surface the inversion symmetry is al- 
ways broken, and the SHG always emerges. It is difficult 
to separate bulk and surface contributions in reflection 
measurement. 

This SHG can arise from various multipole contribu- 
tions, among which an electric-dipole one is the most 
dominant in general. This arises from the third nonlin- 
ear dielectric constant xtjk defined as 

P,{2u) = x%k{'^uj-uj,uj)E,{u)Ek{u;) (27) 

and is nonzero only if the spatial inversion / is broken. 
Thus this electric-dipole SHG appears either at a bulk 
with broken inversion symmetry or at surfaces. Depend- 
ing on magnetic symmetry of the compound, the SHG 
emerges in a different way by varying the incident di- 
rection or polarization of light, and by changing which 
polarization is measured. Therefore, by measuring the 
SHG intensity for various experimental settings, mag- 
netic symmetry of the material can be almost uniquely 
deter.mmpj±^s the recent examples of YMnOa and Cr203 

showNBElEl. 

CaBg is not an exception. Two candidates for its mag- 
netic symmetry, (C-1) 4/m'mm and (D-1) 4'/m'm'm, 
can be distinguished by the SHG, as we shall see be- 
low. Because CaBg breaks inversion symmetry according 
to our prediction, this electric-dipole SHG is nonvanish- 
ing in the bulk. The third nonlinear dielectric constant 
xtjk is decomposed into two parts xtjk = X-jfc + x'jfc\ 
where Xijfc Xijfe' ^-re i-tensors (invariant under time- 
reversal R) and c-tensors (change sign under R), respec- 
tively. In CaBg, the _R/-invariance yields Xijfc — 0- 
The magnetic point-group symmetry of the system de- 
termines which components of the tensor x'ljk van- 
ishing. We can calculate the intensity of the SHG for any 
choice of incident direction and polarization. Among the 
various choices the simplest one is the light propagating 
along the z-axis, i.e. the c-axis in the tetragonal crys- 
tal. As shown in Appendix however, the SHG cannot 
arise in this setting. Thus, the samnle must be tilted 
to observe the SHG, as in YMnOsHO. The simplest 
choice of tilting of the sample is around the [100] axis. It 
is, however, not convenient for the identification of mag- 
netic symmetry for CaBg ; the difference between the two 
possibilities (C-1) and (D-1) is not conspicuous with this 
choice. Instead, we shall tilt the crystal around the [110] 
axis while the incident dicection is fixed as A;||z. Let 
(j) denote the tilting angleEj, which is fixed throughout 
the whole measurement. The details of the calculation 
is presented in Appendix |^. Among various choices of 
polarization of the incoming light E(uj) and that mea- 
sured by the detector, we find the following; when we 
wish to see the difference between the two cases (C-1) 
(D-1), the best way is to measure the SHG intensity for 
the polarization perpendicular to that of the incoming 
light E{uj). Let P{2uj)± denote the component of P{2uj) 
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perpendicular to the incident polarization. Then we get 

eoSin0 (^^e(c)T^2-2 



(C-1) ■.P^=E, 



El + El 



-XliiiEl + El cos^ 4>) - xliiEl sin2 cfj , (28) 
'"''"^ %l'^AEl~ El) cos<^ 



(D-1) -.P^^Ey 



El + El 



-2xliiEl cos </> - xlil E,,Ey sin^ </>) . (29) 

The SHG intensity is proportional to Pj. Thus we pro- 
pose the following experiment for distinguishing (C-1) 
and (D-1). As the polarization of the incoming light is 
rotated, and polarization measured by the detector is ro- 
tated accordingly as is perpendicular to E{uj), we should 
measure the change of the SHG intensity. In (C-1) the 
SHG intensity vanishes when E^ = and in (D-1) it van- 
ishes when Ey = 0. Thus, in this measurement we can 
distinguish between these two. This SHG is observed 
only below T^, which will be an evidence of magnetic 
ordering in the parent compound CaBg. 

Let us consider an effect of domain formation on the 
SHG. For the (a)-type domains, x's change sign for dif- 
ferent domains; this does not, however, affect the SHG 
intensity, which is proportional to . Thus the type (a) 
domains will do no harm. On the other hand, the (b)- 
type domains will kill the difference between (C-1) and 
(D-1), as explained below. If the tetragonal distortion is 
along the x- or y-axis instead of the z-axis, we can calcu- 
late in the similar way as above. The resulting form of P± 
is so complicated that we do not reproduce here. Unlike 
( p8|) (|29|), P± cannot be factorized, i.e. there is no node 
when the incident polarization is rotated. If this type 
of domains are mixed along the direction of light prop- 
agation, the SHG intensity does not have no node, and 
we can no longer distinguish between (C-1) and (D-1). 
Therefore, we should make sure that this (b)-type do- 
mains are not mixed. This can be verified by an untilted 
{(j) ~ 0) setting. If there is only a single domain, we can 
use (^8|)(|2^) and the SHG intensity vanishes when ip — 0. 
An existence of the (b)-type domains will give nonzero 
intensity. To avoid the mixture of these (b)-type do- 
mains, a choice of substrate material will be important. 
It is expected that a good choice of substrate material 
for thin-film fabrication will uniquely set the direction of 
tetragonal distortion to be perpendicular to the thin-film. 

In cases such as Cr203, contribution from magnetic 
dipole cannot be neglected. Let us, therefore, consider 
this contribution for CaBg. This comes from the SH 
magnetic-dipole tensor x™^ defined by 

M,(2l^) = '-^x7;ki2io;Lu,Lo)E^{cu)Ek{Lo). (30) 
n 

Since this is an axial tensor, the i?/-symmetry in CaBg 
allows only the i-tensor component xJ^h^- Hence it per- 
sists both below and above the Neel temperature. Nev- 
ertheless, we show in Appendix |^ that when we measure 



the polarization of the SHG perpendicular to E{uj), this 
magnetic-dipole contribution never appears. This is an- 
other merit of our experimental setting. ___ 

The AF domain topography by the SHGoLSeJ is also 
possible in this material. Domain topography utilizes an 
interference between the c-tensor x*^*^) and the i-tensor 
component x'*''- Thus we should have the SHG from 
both contributions at the same time for domain topogra- 
phy. At least the detected polarization should not be per- 
pendicular to the incident polarization. To obtain better 
contrast for two types of domains, the polarization of 
the irpidcnt light and that of detected light should be 
tunedEil. 



V. DISCUSSION AND CONCLUSION 

Now we discuss the relevance of our results to var- 
ious experiments on CaBg and Cai-ajLaajBg. We be- 
lieve many of the novel magnetic properties can be in- 
terpreted as the ME effect. The ferromagnetism in the 
thin-film CaBg is interpreted as caused by an electric 
field between vacuum and the substrate. For the powder 
sample experiment and the La-doping experiment, the 
explanation is more delicate. We believe that the car- 
riers by La-doping is trapped by impurities/defects and 
create local electric fields. Therefore in these cases, an 
internal electric field and/or a gradient of a strain has 
a random direction, and hence the magnetic moment is 
induced locally due to this mechanism. Without an ex- 
ternal magnetic field, they almost cancel with each other, 
giving zero or quite small uniform magnetization, which 
appears to contradict with the experiments. In fact it is 
consistent with experiments. The compound shows hys- 
teresis, which is thought as an experimental evidence for 
"ferromagnetism". This is quite opposite to our intu- 
ition. The present compound is AF, which usually does 
not show hysteresis behavior in the M-H curve; it is a 
novel feature of doped AF magnetoelectrics. Hysteresis 
implies that the free energy has a double minimum as a 
function of magnetization. In a first sight it is unique to 
ferromagnets; we propose that this is also the case for 
Cai_a;Laa;B6, which we claim is not ferromagnet in the 
bulk. In this material, there is a local electric field in- 
duced by doped impurities. When we vary the magnetic 
field, the two AF domains will switch to each other to 
minimize the free energy. As we explain below, Gibbs 
free energy acquires two minima as a function of M in 
the presence of local electric field, which causes hystere- 
sis. Considering that the ME tensor G is proportional 
to the exciton order parameters Imrj, we can write the 
Gibbs free energy F as 



F = -C{linr])PM 



1 



1 



2x 2a 
B{lmr] f - DPRei], 



p2 - A(lTJir]f 

(31) 



where A, B, C, x, ct are positive constants and D is a con- 
stant. Here, for simplicity, we omitted subscripts for P, 
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M, 7] without loss of generality. By minimizing F in 
terms of Im-q, we obtain Imr; ~ sgn(P A/) A/2B, pro- 
vided C is small. When we substitute it into F, we get 



DPRer]. (32) 



Near the doped impurities, the polarization P is induced 
by a local electric field. In this case, P, M, and 77 has spa- 
tial dependence and should be written as P{x), M(a;), 
T]{x). Although the sign of P{x) changes spatially, one 
can basically regard |-P(£c)| as a constant. In the pres- 
ence of P, this Gibbs free energy has two minima as a 
function of M, implying hysteresis. Note that this is not 
a genuine ferromagnet; a uniform M{x) does not emerge 
in the absence of an external magnetic field, because the 
uniform M{x) costs an elastic energy (VIm77(a;))^. On 
the other hand. If there is no impuritiy, i.e. P — 0, F 
does not have double-minimum structure, and no hys- 
teresis results. Hence, hysteresis in Cai-^LaxBe can re- 
sult from the ME effect together with the local electric 
field near doped impurities. 

Because Cai-ajLa^^Bg is not insulating, the electric 
field is screened. The Thomas-Fermi screening length 
is estimated as 




lOA, 



(33) 



where uq is the Bohr radius_and we used the values e ~ 6 
for the dielectric constaniEj and no ^ (1-4) x 10^'*m~^ 
for a density of electronst3. The appearance of magnetic 
moments via the ME effect is therefore localized within 
this screening length of the impurities or surfaces. 

Other peculiarities of Cai_a;La2:B6 can also be ex- 
plained as well. The high Curie temperature 600K) 
is nothing but a Neel temerature of the parent com- 
pound CaBg, and is not contradictory with a tiny mag- 
netic moment. A rather narrow range (x ^ 0.01) of La- 
doping allowing ferromagnetism is attributed to fragility 
of excitonic order by a small amount of impuritiesLllu. 
Moreover, our scenario is also consistent with, the ex- 
perimental results that deficiency, in Ca-siteso or dop- 
ing of divalent elements like Balil or SiO induces fer- 
romagnetism. It is also confirmed numerically by a 
supercell approach that, imperfections and surfaces can 
induce local momentsE^I. It is hard to explain them 
within the spin-doping scenarioj. Furthermore, strangely 
enough, it is experimentally hard to find a correlation 
between magnetism and electrical resistivity, as soen 
in magnetizaionE2l and in nuclear magnetic resonanceE3. 
This novelty can be considered as natural consequence of 
our scenario; electrical resistivity should be mainly due 
to doped carriers by La, while the magnetization is due 
to local lattice distortion and/or electric field. 

The ESR experiments by KuniiEj also support the 
above scenario. The ESR data show that in a disk-shaped 
Cai-ajLa^Bg {x = 0.005), the magnetic moment only ex- 
ists within the surface layer. Furthermore, the moment 



M does not orient in the direction of H, i.e. it feels 
strong magnetic anisotropy to keep the moment within 
the disk plane. This might be due to the long-range dipo- 
lar energy, and not due to the above scenario. Never- 
theless, it is unlikely that the long-range dipolar energy 
causes such a strong anisotropy. This point requires fur- 
ther experimental and theoretical investigation. Let us, 
for the moment, assume that this strong anisotropy is 
mainly caused by the exciton condensation and the ME 
mechanism. Since this electric field should be perpendic- 
ular to the plane, the strong easy-plane anisotropy paral- 
lel to the surface implies that MJ-E. Therefore, among 
the four cases, (C-1) or (C-2) are compatible, i.e. the 
excitons corresponding to J7(ri5, 1) and J7(ri5, 2) simul- 
taneously condense. Considering the Raman scattering 
data, we conclude that (C-1) is the only possibility com- 
patible with the experiments. 

We briefly mention the relationship between our theory 
and the theories in Refs. ||,^. In the absence of doping, 
our theory is consistent with Refs. |^,^. The difference 
lies in the mechanism of ferromagnetism in doping. We 
proposed that defects, impurities or surfaces induce fer- 
romagnetism. This is different from Refs. |^,^, in which 
ferromagnetism is due to spin alignment of doped car- 
riers. We should note that these two mechanisms can 
coexist. We can distinguish these two contributions for 
ferromagnetism by systematically changing the valence 
of doped impurities. In doping with divalent doping, no 
carriers are doped, and the moment is due to our sce- 
nario. Difference between trivalent- and divalent-doping 
corresponds to the scenarios in Refs. ||,^ Our theory is 
treating the dilute doping limit. We also note that be- 
cause our theory is based on the GL theory, it is treating 
an instability in slight doping. Thus our conclusions do 
not hinder an appearance of rich phase diagrams, as pre- 
sented, for example, in Ref. |^,[^, p"7[p^ , p^ . 

We mention here a role of the spin-orbit coupling. In 
the absence of the spin-orbit coupling, the GL free energy 
(E"*^^) for the imaginary parts of the order parameters is 
written as 

$(^)=A^"^(Im,y,(ri5,+))^ 

i 

+ ^(im ^,(ri5, -)r + >^rL E(i"^ ^.(^25))^ 



^ > aL"' = Af;^ from 



from ( p2| ) in Appendix ^. As Ajj 

the aforementioned microscopic calculation, the conden- 
sation of excitons in Fis and those in occur simulta- 
neously. Furthermore, inspection of ^'■-^"^ in the absence 
of the spin-orbit coupling shows that the quartic-order 
terms in the GL free energy do not lift this degeneracy. 
Sixth-order terms will lift it, and a resulting state will 
belong to either 4/m'm'm' or 4'/m'm'm. The details are 
presented in Appendix Both of them still lead to 
the ME effect in the absence of the spin-orbit coupling; 
this ME effect must be generated from an orbital motion. 
Thus, the AF state in CaBg has an orbital nature as well 
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as a spin nature. With spin-orbit interaction, these two 
are inseparably mixed together. l— , 

Recently, similar noveLierromagnets such as CaB2C2E2l 
and rhombohedral CgoEj have been discovered. They 
share several properties with CaBg, i.e. high Curie tem- 
perature, smallness of the moment and lack of partially- 
filled d- or /-bands. They might be explained by the sim- 
ilar scenario as in CaBg, and indeed one of the authors 
explained the novel ferromagnetisrn—in CaB2C2 within 
the scenario of exciton condensation^. 

In conclusion, we have studied the symmetry proper- 
ties of the excitonic state in the parent compound CaBg, 
and found that the triplet excitonic state with broken 
time-reversal and inversion symmetries offers a natural 
explanation, in terms of the ME effect, for the novel fer- 
romagnetism emerging in La-doping or thin-film fabrica- 
tion. This scenario can be tested experimentally by mea- 
surements of the ME effect and the optical non-reciprocal 
effect in single crystal of the parent compound CaBg. 
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APPENDIX A: DERIVATION OF THE 
REFRACTIVE INDEX - OPTICAL 
NONRECIPROCAL EFFECTS - 

Nonvanishing components for each tensor of the opti- 
cal nonreciprocal (NR) cflfects are determined from each 
magnetic point group asE3 

(C-1) 



(D-1) 



^xx - 
Gxy - 
^xxz 
^zxx 
^xx - 
Gxy ' 
^xxz 



^yy: 

Gyx , 
^XZX 



yzy^ 



zzz ' 
Oizz, 



yyz 



yzy^ 



In calculating the optical properties in transmission or in 
refiection, the tensor Aaf3-y defined as 



= —tfi-isGaS — ^ajsGps + -Lo{a'^p^ -I- a^„^)(Al) 



plays a central rolecfl. Its nonvanishing components 
for each magnetic group are summarized in Table 2 of 
Ref. ||. They are 



-Gxy + -^^^i^'zxx + ^'xxz)^ 



(D-l) 



Ax 

Az 



Ayyz 2Ca: 



y 



OJCC 

-A — — 4 — ~ A 

' -^xzx -^zyy ^ 



1 



'Gxy + -^i^ia'zxx + '^'xxz)y 



AxXZ Ayy2^ IGxy 



The refractive index n for each polarization is obtained 
as a solution of the following equation. 



n'^{aa(Jl3 - Sap) + Safj + Eg + C^o"-Cr^^Q/37 



■En 



0. 



(A2) 



For incident direction along the ±z-direction cr — k/k — 
(0,0,0-2) — (0,0, ±1), an equation determining the re- 
fractive index n is 



+ Cq + c^onC5(7z \ / Ex 

(C-1): I -n^ + e-^x + c^l„nC5CT, ^ \{Ey\^0, (A3) 

e^^ez + ciiquCqcFz J \ 

-n^ + Cq^Ex + cfionC^a^ \ / Ex 

(D-l): I -n^ + e^hx - cfionC^a, \ \ Ey \ = 0, (A4) 

E._ 
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eo + Oizz C5 



where = eo + "^x = £0 + ^1 
2Gxy + Lua'^^^, and Cg = ^o-'zzz- 
refractive index n is given by 



(C-1) : [E\\x, E\\y] ■ n ^ ^ [ (Jz^Cs + \/ f^lCl + 4/ioe. 



By using thjem, we can derive the equations for the SHG, 



Thus for both cases the as in Cr203c2l. The nonhnear polarization P(2lS) induced 

by X%k is 



(D-1) : {E\\x\ 



- [ azfJ-oCs 



[E\\y] ■ 2^ -CTz/^oCs 



for cr ^ (0,0, ±1). 

In order to consider an effect of domain formation of 
the (b)-type, i.e. domains with different directions of 
distortion, let us study the light propagating along the 
x-direction in the similar way as above. An equation for 
(T=(a,,0,0) = (±1,0,0) is 



6n 



^0 




for both (C-1) and (D-1), where we set C7 = 
^uj{a'^.j.x -f a'xxz)- Thus the refractive index is 



-Gx 









f 


'^'XxxzExEz ^ 




(C-1) 


P{2uj) 


= eo 




2x"MEyEz (] 










\ Xzxx 


iEl + E^y) + xtilE^J 










} 


'^Xxxz^x-^z ^ 




(D-1) 


: P{2lo) 






-2xli'lEyEz (] 


52) 








e(c) 
\ Xzxx 


iE^x-E^y) + xtilE^J 





These contribute to the source term S"^ for the SHClEl as 



d^P{2uj) 



-Auj^HaP{2uj)- 



[E\\v] : n 



£0 £2:, 



\E\\z\ : n=xle-\z ( 1 + 



for (T = (±1,0,0). The latter one with E^z contains a 
longitudinal component of E, i.e. parallel to a;. It is 
called an S-wave (skew wave) in Ref. and can only 
propagate only inside the sample. For cr — (0, ±1, 0), we 

get 



[E\\x] : n 



\E\\z] : n=,/e^hz ( 1 + ^Cg 



These imply that for (t\\x and cr ||y the NR optical effects 
do not take place. 



Therefore, light propagating along the z-axis, i.e. the 
c-axis in the tetragonal crystal, S is along the z-axis 
and the SHG cannot be generated. Thus, the inci- 
dent disectipn must be tilted to observe the SHG, as in 
YMnOaEj'H. For practical calculations, it is more conve- 
nient to tilt the crystal while fixing the light propagation 
along the z-axis. This kind of treatment is also adopted 
in Ref. |l|. The simplest choice of tilting of the sample is 
around the [100] axis, but this choice does not manifest 
clearly the difference between (C-1) and (D-1). Instead, 
we shall tilt the crystal around the [110] axis while the 
incident direction is fixed as fc||z. We fix the tilting angle 
(j), and we fix its value throughout the whole measure- 
ment. The calculation is similar to the one for YMnOa 
in Ref. |5l|. The new coordinate system (A'^, A'^, A'^) fixed 
to the crystal is related to the original one {A^, Ay, Az) 
as follows; 



APPENDIX B: CALCULATION OF THE SHG 
INTENSITY TO DISTINGUISH (C-1) AND (D-1) 

The nonvanishing components of the nonlinear 
electric-dipole tensor Xij'k (polar c-tensor) are given asEl 



(C-1) 


A.ZZZ 


/Vzxx 


Azyy I 






A.XZX 


- y"(") 

A.XXZ 


- y"(") - 

Ayzy 


Ayyz I 


(D-1) 


ye(c) 

A.ZZZ 


y^zxx 


- ^yeic) 
Azyy 






yeic) 

A.XZX 


A.XXZ 


- ^yeic) 
Ayzy 


Ayyz 



A'x 

A'y 1 = 

A' 



1 

V2 




First, the polarization E{ijj) = (E ^ , Ey , 0) is transformed 
to the new coordinate. Then (Bl) (B2) are applied to get 
the nonlinear polarization P(2a;), and transform it back 
to the original coordinate {Ax,Ay,Az). The resulting 
form for the nonlinear polarization is 
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(C-1) : P{2uj) = £0 



(D-1) : P{2u;) = £„ 



-^Xx'^xlExEy sin I 



COS © sm 



- Xzx 



■^^\El+ El COS' 



■x\^'. 



- C0S2 



-2x'xxIeI coS(?!)sin( 



0) sin (j) - x'zzIeI s\v? 
cos (j) + Xzi'^zEy cos sin^ 



X^zxl)E^Ey cos (/» sin (/> - XzHe^ sin^ 

: J. I o,e(c)j;i 7;n 2 J, I ,e(c)77.2„- 



-2xxxlExEy SIT? (j) + 2x''zxlE^Ey cos^ (/) + x'z^zIe'I s\v? (j) cos ( 



(B4) 



(B5) 



In order to see the difference between (C-1) and (D-1), 
we find that the best way is to measure the SHG intensity 
for the polarization perpendicular to that of the incoming 
light E{uj). The projection of P{2uj) onto the direction 
perpendicular to both E{u}) and 2-axis is just P(2w)i = 
P{2uj) ■ n±, where n± ^ E~'^{-Ey, E^, 0), Eq = 



El 



(C-1) 



(D-1) 



£2. Thus 



Pi = 



El + El 



i£l{El+ El cos 



' X sin {2x1^^1 El sin^cj, 



P± = 



l^^^ El sin' A ,(B6) 



■ySincf>{2x%ii{El-E',)cos^ 
^xliiEl cos 4, - xfzlExEy sin' A . (B7) 



El + El 



These contribute to the source term 5™ for the SHOE 
as 



dM(2uj) 
dt 



- Ey sin (/) cos Ey 



8w 
6 



,2 xl'^^ Ey sin (f) COS (j)E. 



Therefore, this S'™ is parallel to the incident polarization 
E{uj). Thus, when we measure the polarization of the 
SHG perpendicular to E{uj), this magnetic-dipole con- 
tribution never appears. 



APPENDIX D: LIMIT OF ZERO SPIN-ORBIT 
COUPLING 



which are identical with (P8| )(p9[). Note that (B6) is pro- 
portional to Ex, while (|B7| ) is proportional to Ey] this 
enables an identification of the true ground state. 



APPENDIX C: MAGNETIC DIPOLE 
CONTRIBUTION TO THE SHG 



Its nonvanishing components are given as 



(C-1)(D-1) 

■ Xxyz Xxzy Xyxz 



\yzx 7 



_ m(i) 
\zyx ' 



By using them, we can derive the equations for the SHG, 
as in Cr203Ej. The nonlinear magnetization M{2u}) in- 
duced by xTjk is 



M{2uj) 



2eoc 



™^5£;„sin 



~Xxy 




In the absence of the spin-orbit couphng, the GL free 
energy should be invariant under the spin rotation; by 
lengthy calculations this leads to relations 

A2 + Ai = 2Ai, ^2 - ^4 = A3/%/2, 

B2 + S4 = 2Bi, B2- B4 = B3/V2. (Dl) 

When we substitute them into $(2) we get 

^(2) ^ ^(Im ry,(ri5, +))' + A^^ ^(Im m{ri5,~)f 

i i 

3 

+ AiJ^ ^(Re 77»(ri5, +))' + A«<= 5](Re r^^iV,,, -))' 

i i 

+A^^:.E(i^e^^(r25))^ (D2) 

where T7(ri5,±) are new bas is functions, defined in ( pT| ) 
(|l|). The coefficients in (|D2|) are 



\Im _ 352 -B4 ,Ini _ \Im _ n 

^+ - 2 ' ~ 1^25 ~ -^4, 

xRo 3^2 - A4 

A+ = 5 , A_ =Ar,5-A4. 



(D3) 
(D4) 
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Eqn.(D2) is invariant under the cubic grou p op erations 
in the orbital space. Thus, each term in ( |D2|) can be 
classified into irreducible representation of the Oh in the 
orbital space, which facilitates subsequent discussions for 
magnetic properties. Focusing on the imaginary parts of 
the order parameters, we get the result 



$(2) = J2 



l—u,v 



where 



(D5) 



transform according to the T[ and r'12 representations of 
Oh, respectively, under the rotation of the orbital space. 
As X^_^ > a'™, the condensation of excitons in F'j^j occur, 
and we shall focus only on r'12 There are two types of 
degeneracies in (D5). One is in the summation over i, 
which reflects the SU(2) symmetry of the spin space. It 
will never be lifted in the absence of the spin-orbit cou- 
pling. The other is in the summation over I, i.e. in the 
direction in the rju-Tjy plane. This is lifted in the sixth 
order in rji. To see this, let us write down higher or- 
der terms. They are written conveniently in terms of a 
complex order parameter w = i]^ + irj^ as 



<^ = K^y\w{T'^2,i)? + K2 



+ K, (^\w{T[2,i)\'^ +K,Re (^[wir\2,i)?^ 
where Kj are constants. This is minimized when 



(I) argw = f n if K4 <0 

(II) argw = f (2n -f 1) if > 0, 

where n is an integer. The order parameters are 

(I) (Im?7,(Q,), Imij.iQy), Imij.iQ,)) = c(l, 1, ^2)P) 
(II) (Im77,(g,), luiTj.iQy), Im?7,(Q,)) - c(l, -1,((I)7) 

where we write down only one among three equivalent 
directions for each case. Both (I) and (II) have tetragonal 
distortion. The magnetic point group for each case is 
(I) 47m'm'm = D4 x {E, RI} (II) 4/m'm'm' = Dzd X 
{E, RJ}- The ME property tensors are easily written 
downEj; 



(I) 



(D8) 



Gi \ /Gi 

Gi , (II) -Gi 
00 Ga/ \0 00 

in the cubic coordinates. Note that the ME effect in the 
absence of the spin-orbit coupling cannot originate from 
spins. This shows that the orbital moment exists in this 
system. 

It is helpful to compare the states (C-l)-(D-2) with 
(I) (II). The forms of the order parameters are given 
in Eqs.(||) (H) for (C-1) and (D-1), and in Eqs.(0) 
(p^) for (I)(II). Naively one may expect that (C-l)-(D- 
2) are included in (I) or (II), because (C-l)-(D-2) assumes 
nonzero spin-orbit coupling, whereas (I) (II) assumes an 
absence of the spin-orbit coupling. It is, however, not 
the case, as seen from the different ME tensors in the 
two cases. The reason is the following. In the presence 
of the spin-orbit coupling, the degeneracy in the GL free 
energy $ is lifted in the quartic order, resulting in (C-1)- 
(D-2). In contrast, without the spin-orbit coupling, the 
degeneracy is lifted in the sixth order, leading to (I) (II). 
Because of this difference in the lifting of degeneracy, the 
realized states are different between the cases with and 
without the spin-orbit coupling. 
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